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ABSTRACT 

Using hydrodynamical simulations of the ACDM cosmology, that include both radia- 
tive cooling and a phenomenological model for star formation and supernovae feedback, 
we investigate the impact of galaxy formation on the X-ray properties of groups at red- 
shift zero. Motivated by the observed "break" in the L^ — T^ relation at kT^ ~ 1 — 2keV, 
our feedback model is based on the assumption that supernovae imprint a temperature 
scale on the hot gas, with the star formation rate and corresponding reheated gas mass 
then depending only on the available energy budget. We demonstrate that a strong 
feedback model with a heating temperature comparable to this break (fcTgN = 2keV), 
and an energy budget twice that available from supernovae (e = 2), raises the core 
entropy of groups sufficiently to produce an adequate match to their observed X-ray 
properties. A lower value of e increases the star formation rate without significantly 
affecting the X-ray properties of groups, and a model with e ~ 0.1 reproduces the 
observed fraction of baryons in stars. However, a heating temperature which is lower 
than the virial temperatures of the groups leads to an excess of cooling gas that boosts 
their X-ray luminosities, due to the failure of the reheated material to escape from 
the gravitational potential. A limited study of numerical resolution effects reveals that 
the temperature of poorly-resolved objects is underestimated, therefore (in our case) 
a fully-resolved group population would lead to a steeper — Tx relation, bringing 
our results into even better agreement with the observations. 

Key words: hydrodynamics - methods: numerical - galaxies: formation - X-rays: 
galaxies: clusters 



1 INTRODUCTION 

Groups ^ are the most common agglomerates of galaxies 
in the low-redshift universe and also harbour the major- 
ity of currently detected baryons (Fukugita, Hogan & Pee- 
bles 1998), observed because of their emission of soft X- 
rays {kT ^ 0.1-1 keV). The simplest explanation for the 
properties of this gas is that it was shock-heated during 
the gravitational collapse and virialization of the system. 
For self-similar (S7 — 1) systems, this leads to a simple 
scaling between X-ray luminosity and temperature, Lx cx 
T^/^(l -I- zf^^A{T), where A(r) denotes the cooling rate of 
the gas. 

In the largest clusters, the virial temperature is 
sufficiently high that the gas predominantly emits 
bremsstrahlung radiation, leading to the scaling Lx tx 
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^ We use the term "groups" to mean systems with X-ray tem- 
peratures of order 1 keV or lower 



(Kaiser 1986). This result was recently vindicated by Allen, 
Schmidt & Fabian (2002) using high-quality Chandra data. 
In smaller systems, the increasing importance of line emis- 
sion boosts the luminosity and so should flatten the Lx — Tx 
relation. Observationally the converse is true, with the slope 
of the relation for low-mass groups being as high as ~ 6 (e.g. 
see Xue & Wu 2000 for a recent compilation). This result is 
also reflected in the low amplitude of the soft X-ray back- 
ground, around an order of magnitude lower than predicted 
from gravitational heating models (e.g. Pen 1999; Bryan & 
Voit 2001; Wu, Fabian & Nulsen 2001). 



Evrard & Henry (1991; see also Kaiser 1991; Bower 
1997; Voit et al. 2002) argued that the deficit in luminosity 
is due to additional physical processes raising the entropy 
of the gas. (Entropy is invariant under adiabatic processes, 
making it a useful variable for characterizing the effects of 
cooling and heating on the gas distribution.) Gravitational 
heating leads to the mean entropy of a self-similar system 
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scaling as s oc T/{1 + z)^ ^. Increasing the entropy of the 
gas corresponds to a hotter but less-concentrated distribu- 
tion. Consequently, for a given increase in entropy, the X- 
ray emissivity is lowered preferentially in low-mass systems, 
steepening the L^ — T^ relation. Ponman, Cannon & Navarro 
(1999) measured the entropy of groups and clusters at 10 per 
cent of their virial radii, and showed that low-mass clusters 
and groups deviated from the self-similar prediction, tend- 
ing to a constant value of ~ lOO/i"^''^ keV cm^ (albeit with 
very large scatter). 

The precise origin of the excess entropy is still unclear, 
but is most probably a direct consequence of the galaxy for- 
mation process (Evrard & Henry 1991). For example, galac- 
tic winds from supernovae explosions and/or active galactic 
nuclei (AGN) could heat the gas, raise its entropy, and allow 
it to convect out to lower density environments. Winds have 
already been observed in low-redshift starburst (e.g. Heck- 
man 2000) and high-redshift Lyman Break (e.g. Pettini et 
al. 2001) galaxies, and the presence of metal-enriched gas 
in low-density regions suggests galactic material has been 
expelled from haloes (e.g. Schaye et al. 2000). Note that 
in hierarchical models, at least some heating (feedback) is 
required in order to regulate the star formation rate (e.g. 
White & Rees 1978; Balogh et al. 2001). 

The amount of energy required to generate excess en- 
tropy depends on environment. If the gas was close to 
mean density when heated (the pre-heating scenario) only 
~ O.SkeV of energy per particle is required (Lloyd-Davies, 
Ponman & Cannon 2000). More realistically, if the heating 
occurred in high-density regions, as much as 1-3 keV per par- 
ticle is required (e.g. Wu, Fabian & Nulsen 2000; Kravtsov 
& Yepes 2000; Bower et al. 2001; Borgani et al. 2002). Note, 
however that the total energy production is smaller in the 
latter case as only gas in the cores of haloes need be heated. 

Excess entropy is also generated by the process of 
star /black- hole formation itself, removing cold, low-entropy 
material and causing higher-entropy material to flow in to 
replace it (Knight & Ponman 1997). Models using only ra- 
diative cooling to generate the excess entropy have been im- 
plemented by a number of authors, and have proved success- 
ful in reproducing a number of X-ray scaling relations (e.g. 
Pearce et al. 2000; Bryan 2000; Voit & Bryan 2001; Muan- 
wong et al. 2001,2002; Wu and Xue 2002a,2002b; Dave et al. 
2002). Voit & Bryan (2001) demonstrated that the entropy- 
temperature locus of gas with a range of temperatures and 
densities set by equating their cooling times to the age of the 
universe, closely matches the relation observed by Ponman 
et al. (1999). Voit & Bryan further argued for the presence 
of an entropy floor by asserting that cooled gas will either 
form stars or be heated (by stars or AGN) and move to a 
higher adiabat. 

The aim of this paper is to investigate how galaxies 
which form in a cosmological simulation impact upon the 
X-ray properties of groups. This treatment allows us to in- 
corporate cooling and heating processes in a self-consistent 
manner. Specifically, we include the effects of both radiative 
cooling/star formation and feedback of energy from super- 
novae in our simulations, but neglect an AGN component 



^ We use the term entropy to mean s = T j p"' ^, and assume 
7 = 5/3, for a monatomic ideal gas 



(although in broad terms, the effects of including AGN may 
largely mimic the desired effects of supernovae, since both 
are ultimately powered by cooling gas). Motivated by the 
observed "break" in the Lx — Tx relation at the interface 
between groups and clusters {kT ~ 1 — 2keV), our feedback 
model is fundamentally based on the assumption that cold 
gas is reheated to a fixed temperature and the mass ratio 
of reheated gas to stars is determined by the overall energy 
budget. As we shall see, our model is able to reproduce the 
required entropy excess, provided that the heating temper- 
ature is sufficiently high to allow the majority of reheated 
gas to escape from haloes, while varying the energy budget 
primarily controls the fraction of baryons in stars. 

We organize the remainder of this paper as follows. In 
Section |5| we outline our method, giving details of the sim- 
ulations as well as our approach to cooling, star formation 
and feedback. Section |3 presents the predicted X-ray group 
scaling relations for our fiducial simulation at z = 0. In Sec- 
tion 0] we investigate our feedback model in more detail, 
concentrating on the effects of varying parameters. In Sec- 
tion |K1 we perform a limited study of numerical resolution 
effects and summarize our conclusions in Section |S] 



2 SIMULATION DETAILS 
2.1 Overview 

Results are presented for the currently-favoured ACDM cos- 
mology, setting the density parameter, f2m = 0.35, Hubble 
constant, h — _ffo/100kms~^ Mpc~^ — 0.71, cosmological 
constant, JIa = A/3-ffo — 0.65, baryon density, 57;, — 0.038, 
and power spectrum normalisation, as = 0.7. The latter 
value, recently measured from both cluster abundances (e.g. 
Seljak 2002) and galaxy clustering (Lahav et al. 2002), 
is around 20 per cent lower than previous determinations 
which deduced erg ~ 0.9 (although, weak lensing statistics 
still prefer the higher value; see e.g. Refregier et al. 2002). We 
find that increasing erg from 0.7 to 0.9 has only a small effect 
on the results presented in this paper, but the lower value 
gives us a numerical advantage, requiring fewer timesteps to 
evolve to 2; = 0. 

Simulation data were generated using a parallel 
(OpenMP) version of the hydra code (Couchman, Thomas 
& Pearce 1995; Theuns et al. 1998), a combination of the 
AP^M algorithm (Couchman 1991) to compute gravitational 
forces and Smoothed Particle Hydrodynamics (SPH; see 
Monaghan 1992) to calculate gas forces. For the gas, neigh- 
bour lists were computed using the recursive binning tech- 
nique described in Theuns et al. (1998), allowing the den- 
sity to be estimated for gas particles in low-density regions. 
SPH forces were calculated using an artificial viscosity based 
on pairwise particle separations (Monaghan 1992) rather 
than the local velocity divergence. Our SPH implementa- 
tion now also includes the force-correction terms suggested 
by Springel & Hernquist (2002), resulting in improved en- 
tropy conservation. 

Initial conditions were generated using the COSMIC 
package supplied with the publicly available version of hy- 
dra. We made a small modification to cosmic so that the 
large-scale structure is preserved when increasing the numer- 
ical resolution. For our main results, we adopted a comoving 



The impact of galaxy formation on X-ray groups 3 



box-size of 50/i~^Mpc, populated with a perturbed grid of 
2,097,152 (128^) dark matter and gas particles each. (Note 
that the box-size is big enough to give statistically meaning- 
ful results; Bryan & Voit 2001.) These choices set our gas 
and dark matter particle masses to rrigas = 6.3 x 10** h~^MQ 
and rridark = 5.2 x 10^ h~^MQ respectively. We have chosen 
this resolution so not to expect significant heating of the 
gas due to two-body encounters with the dark matter par- 
ticles (Steinmetz & White 1997). Gravitational forces were 
softened on short scales using a spline kernel with an equiva- 
lent Plummer softening length of tp = es/2.34 = 10/!.~^kpc. 
Gas forces were also limited by fixing the minimum SPH 
smoothing radius to the gravitational softening 2/imin ~ ^s- 
All but the smallest runs were performed on the UK 
Cosmology Consortium COSMOS Origin 3800 supercom- 
puter in Cambridge, using between 10 and 20 processors. 
Simulations were started at z = 49 and evolved to 2: = 0, 
typically taking between 2000 and 4500 steps, depending on 
the choice of resolution and feedback model. 



2.2 Radiative cooling and star formation 

Radiative cooling was implemented following the method of 

Thomas & Couchman (1992), using tabulated cooling rates 
given by Sutherland & Dopita (1993). For the present paper, 
we adopt a global metallicity scaling law, Z = 0.3{t/to)ZQ, 
to approximate the gradual enrichment of the intergalac- 
tic medium. In principle, we could model metal enrichment 
sclf-consistcntly, since it is a by-product of the supernovac 
explosions. (We have already done this to study the impact 
of dwarf galaxies on the high-redshift intergalactic medium; 
Theuns et al. 2002.) In order to simplify matters however, 
we have decided to leave such a study to a subsequent paper. 

We follow the star formation prescription used by Kay 
et al.(2002, hereafter K2002), converting gas with 5 > 5000 
and T < 50, OOOK into stars at the end of each timestep. 
This density threshold is low compared to measured empir- 
ical density thresholds in galaxies (e.g. Kennicutt 1998) but 
is close to our resolved density maximum. Increasing the 
density threshold above the resolution threshold induces a 
delay to the star formation, increasing the gas to stars ratio 
within each galaxy (K2002). 

Subsequent to the star formation routine, we identify 
groups of 10 stars within the current gravitational softening 
length and merge them into a single galaxy fragment. These 
fragments are comparable in mass to the dark matter par- 
ticles and the reduction in the number of particles provides 
a modest decrease in the CPU time per step. 
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Figure 1. Cumulative number of supernovae per solar mass of 
stars formed versus time since the star formation episode, using 
power-law fits of Lia et al. 2002 and assuming a Salpeter IMF. 

solve explicitly for the flux of material between each phase, 
using a set of physical arguments as laid out by McKee & 
Ostriker (1977). This approach has been implemented and 
proved effective by various authors, using both grid-based 
codes (e.g. Yepes et al. 1997) and SPH (e.g. Hultman & 
Pharasyn 1999; Springel & Hernquist 2003). The main ad- 
vantage of this method is that it does not depend explic- 
itly on numerical resolution, and so should be less sensitive 
to discreteness effects. On the other hand, the phases can- 
not move independently (although see Springel & Hernquist 
2003). 

An alternative approach, which we take here, is to allow 
each particle to represent only one phase and use a set of 
simple rules to overcome any associated discreteness prob- 
lems. This has the advantage of being much simpler to im- 
plement than the previous method and allows each phase to 
move freely. K2002 reviewed a number of feedback models 
within this context; in this paper, we consider an improved 
version of the thermal feedback approach studied by K2002 
and outline details of this method below. 



2.3 Feedback 

The essence of any feedback model is to create a diffuse 
hot phase which powers a wind, providing a mechanism to 
transport gas away from star forming regions. However, ex- 
plicitly modelling a multi-phase medium in a cosmological 
simulation is impossible because of its limited resolution. 
This problem has prompted two approaches to implement- 
ing feedback within such simulations. 

One method is to assume that each resolution element 
(in the case of SPH, a single gas particle) is representative of 
a multi-phase medium (cold and hot gas and stars), and to 



2.3.1 Feedback energetics 

The amount of energy available from each star particle at 
time, t, to heat the gas over a timestep. At, is 

SsN = 10^1 e (^^^ ^ ^ ° [na(r)/2 + rii(r)] dr erg, (1) 

where m* is the mass of a star particle, to is the star forma- 
tion time and ria & rn arc the rates of Typo la and Type II 
supernovae (SNe) per solar mass of stars formed. (We as- 
sume each Type la supernova contributes half the energy of 
a Type II supernova.) 
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Lia et al. (2002) tabulated r{r) for various initial 
mass functions (IMF), using the Padova stellar evolution- 
ary tracks (Portinari et al. 1998). We adopt the Lia et al. 
rates for a Salpeter IMF (given in their Table 2) and illus- 
trate integrated versions in Fig. The contribution from 
Type II SNe rises dramatically after approximately 3 mil- 
lion years, lasting for around 70 million years. Type la SNe 
start contributing later, at around 100 million years, but 
are still going off after 10 billion years. The total number of 
Type II SNe is around an order of magnitude higher than 
the number of Type la SNe, and so provides the dominant 
source of feedback of energy to the interstellar medium. 

The parameter e modulates the overall energy available 
from the SNe, which we term the feedback efficiency. Thorn- 
ton et al. (1998) suggested that the maximal value (e — 1) 
may be at least a factor of 5 too high once radiative losses 
are accounted for. As we shall see in Sectional our model 
does not necessarily require a high value of e to generate the 
required excess entropy; varying this parameter primarily 
alters the fraction of baryons in galaxies. 

2.3.2 Feedback implementation 

The thermal feedback model studied by K2002 is effective in 
reducing the masses/luminosities of galaxies but fails to re- 
produce the X-ray scaling relations of groups. In their model, 
each star particle distributes the energy available at each 
timestep (equation over the surrounding gas. For equal 
mass particles, the maximum amount of energy a gas parti- 
cle can receive from a single star particle is enough to heat 
it by ~ 2keV (for e = 1). In practice, gas particles escape 
before receiving this amount of energy and the result is an 
average heating temperature around an order of magnitude 
lower. Much of this gas then fails to escape from group and 
cluster-sized haloes with comparable or higher virial tem- 
peratures, and instead re-cools, boosting their X-ray lumi- 
nosities (we will demonstrate this effect in Section |1J. 

We have devised an alternative method with the clear 
intention of heating the gas to a higher temperature without 
necessarily increasing the available energy. In the context 
of X-ray scaling relations, our new model is motivated by 
the presence of a "break" in the Lx — Tx relation at fcTx ~ 
1 — 2keV; such a feature may be an imprint of a temperature 
scale introduced by the heating mechanism. We base the 
model on the physical assumption that as more and more 
supernovae go off, it is not the temperature of the gas that 
is progressively increased (as in the previous method), but 
the mass fraction of material converted from the cold phase 
to the hot phase. To implement such a model, we equate the 
available energy from a mass of newly formed stars, m* to 
the energy required to raise the temperature of a mass of 
cold gas, 5mgas by a fixed amount, Tsn 

(5mgaB (3fcrsN/2/xmH) = -Esn, (2) 

where Ssn is given in equation^and we set fimn = 10~^*g, 
the mean molecular weight of a gas with Y — 0.24, Z = 0. In 
practice, Tsn is much greater than the temperature of the 
cold gas (~ lO'^K) so the heating introduces a temperature 
scale to the hot gas. 

Assuming now that m* is the mass of a star particle, 
the fractional mass of a gas particle capable of being heated 
to Tsn is 



using the fact that m* — mgas in our method and R{t — 
to, At) is the integral on the right hand side of equation Q 
(For a typical timestep. At ~ 5 x lO'^yr, (R) ~ 0.001.) Note 
that in K2002's model, /heat ~ 1 and feTsN is calculated, the 
opposite of our model. For every star particle, we calculate 
the radius which encloses 32 gas particles, and (since gener- 
ally /hoat < 1 ) apply equation 13 to the first gas particle 
within this radius with T < 50, OOOK. Since each gas parti- 
cle represents only one phase in our simulations, we can only 
heat the whole particle to a temperature Tsn or not heat it 
at all. To mimic the effect of heating a fraction of the gas, 
we draw a random number for each particle, < r < 1, and 
heat it to Tsn if /heat < r. 

For our fiducial model, for which we present results in 
the next section, we set e = 2 & feTsN ~ 2keV, twice the 
maximum available energy from standard supernovae calcu- 
lations, but in reasonable agreement with that suggested in 
previous studies (e.g. Kravtsov & Yepes 2000; Wu, Fabian 
& Nulsen 2001; Bower et al. 2001; Borgani et al. 2002). As 
we will show, this model predicts X-ray properties of groups 
in good agreement with the data. We will demonstrate in 
Section |1| however, that a good match to the scaling rela- 
tions is reasonably insensitive to this choice of e, but does 
require the heating temperature to be sufficiently high in or- 
der for gas to escape from the gravitational potential of the 
halo. The fiducial temperature, which is significantly larger 
than typically found when using K2002's model, also alle- 
viates the need for preventing gas from cooling radiatively 
over a fixed period of time (as implemented by K2002), be- 
cause cooling gas only starts to become thermally unstable 
at lower temperatures (feT ~ 0.2- IkeV for Z = 0-0.3Zq). 



3 X-RAY PROPERTIES OF GROUPS IN THE 
FIDUCIAL MODEL 

We begin by demonstrating that our fiducial feedback model 
provides an adequate match to current observational data at 
z = 0, given the level of uncertainty in these measurements. 

Our group catalogues are constructed using the proce- 
dure outlined by Muanwong et al. (2002, hereafter M2002), 
where full details may be found. To summarize, spheres of 
baryon and dark matter particles are located with constant 
mean internal overdensity, A, relative to the comoving crit- 
ical density. For our chosen cosmology, A = 111 is approx- 
imately the predicted value by the spherical top-hat col- 
lapse model for a virialized sphere (Eke, Navarro & Frenk 
1998). Unless specified, we set A — 200 for our main results 
but note that the X-ray properties are insensitive to using 
A — 200 rather than A — 111. Only objects containing at 
least 500 baryon plus 500 dark matter particles are retained, 
setting a minimum halo mass of 2.9 x 10^^/i~^Mq. At z — 0, 
our fiducial catalogue contains 172 objects. 

Once the catalogues were generated, it was found that 

^ Exceptions to this rule can occur if e/fcTgN > 5; this is true 
for only one model considered in this paper, where e = 1 & 
^^SN = O.lkeV. As we shall see however, this particular model 
fails because fcTgN is too low, not because less gas is reheated. 



The impact of galaxy formation on X-ray groups 5 







, , , 1 






_ 
- 


















O 






/- 










Or. y/ 








o 


o 


/ 












/ 








% 




^y / 












/ 












/ 










/ 

/ 












/ 

/ 








mi ^ 


























y. o /' 












/ 












/ 












/ 










/ 

/ 

/ 

/ 

/ 

/ 

/ 

, / , 11 , 


/ 

, , , 1 









10.0 



1 .0 

kT./keV 



10.0 



0.1 

0.1 




.0 10.0 100.0 

M„/10"h"'M, 



Figure 2. Core entropy (measured at O.lrvir) against X-ray tem- 
perature. Stars represent groups in the fiducial model and dia- 
monds from the Radiative simulation in M2002. The solid line 
represents the trend observed by Ponman et al. (1999) and the 
dashed line is a fit to the self-similar relation obtained from Non- 
radiative simulations. 



a small fraction (~ 5 per cent) of groups have X-ray tem- 
peratures very much in excess of the average for their mass 
(i.e. there are hot outliers on the temperature-mass rela- 
tion). On inspection, these objects are infalling sub-clumps 
on the outskirts of larger groups, but have hot gas which 
has been shocked to the temperature of the more massive 
object. We choose to eliminate these objects by setting the 
requirement that each group in the catalogue has an X-ray 
temperature less than twice the mean value for its mass. 



3.1 Entropy 

The effect of the removal of hot gas from a halo on its X-ray 
properties can best be characterized by the change in en- 
tropy of the gas. As the hot gas is depleted, higher entropy 
material flows inwards and is adiabatically heated (i.e. it de- 
velops a hotter, less concentrated distribution) in order to 
maintain gravitational support. This leads to a rise in the 
X-ray temperature of the gas but a decrease in its X-ray 
luminosity (Knight & Ponman 1997; Pearce et al. 2000). 
Ponman et al. (1999) demonstrated that lower temperature 
objects have shallower surface brightness profiles and con- 
sequently larger excess entropies, relative to the expected 
self-similar relation for a model neglecting galaxy formation 
(we refer to this as a Non-radiative model). 

Following Ponman et al., we measure the entropy of 
each object at O.lr-vir as a function of X-ray temperature. 
The latter is evaluated within a soft-band of energy range 
0.3-1.5 keV 



J2i mipiKsoft{Ti, Z)Ti 
m,piksott{Ti, Z) 



(4) 



Figure 3. X-ray temperature-mass relation for the fiducial model 
(stars). Temperatures calculated when excluding particles from 
within the cooling radius are shown as squares. Also plotted is 
the M2002 Radiative sample (diamonds). The dotted line is the 
relation between virial temperature and mass. Dashed and dot- 
dashed lines represent best-fit relations to the observational data 
presented in Xu, Jin & Wu (2001) for NEW & /3 model surface 
brightness profiles respectively. 



lO^K) within the virial radius of the group, each with mass, 
rui, density, pi and temperature, Ti. The soft-band cooling 
function, Asoft is calculated using the Raymond & Smith 
(1977) models. We then estimate the entropy of each group 



where 



r) 3M 
TTiH 4:TTr^5r 



(5) 



(6) 



where the sum is performed over all hot gas particles (T > 



is the electron density, averaged over all hot gas particles 
(of total mass, M) within a shell of width, Sr = O.lrvir, 
centered on r = O.lrvir- Only groups containing 10 or more 
particles within this region are plotted. We set the number of 
electrons per baryon, ri — 0.88, appropriate for a primordial 
gas {Y = 0.24, Z — 0) although insensitive to the assumed 
range of Z. 

Fig. 121 shows the simulated entropy-temperature re- 
lation (stars), as well as groups/clusters from the Radia- 
tive simulation studied by M2002 (diamonds). Our fiducial 
model agrees well with M2002's result where the two samples 
overlap, as well as with the general trend observed by Pon- 
man et al. (solid curve) , clearly demonstrating an entropy in 
excess of the self-similar prediction (sx oc Tx, dashed line). 
The increasing impact of the galaxies with decreasing halo 
mass leads to an increase in the excess entropy, consistent 
with the expectation of an entropy floor at Sx ^ 100 ft"'^'^^ 
keV cm^. 
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3.2 X-ray temperature 

Increasing the entropy of the gas leads to a higher X-ray 
temperature at fixed mass, resulting in an increase in the 
normalisation of the temperature-mass relation (e.g. Muan- 
wong et al. 2001; Thomas et al. 2002). We plot the fiducial 
temperature-mass relation in Fig. |3 (stars), with the solid 
line representing a least-squares fit to the data. The square 
points show the result of excluding emission from within the 
cooling radius* of each object. This has the effect of slightly 
increasing the temperatures of some of the larger objects, 
without significantly affecting the slope or normalisation of 
the relation. 

The result is in reasonable agreement with the Radiative 
model studied by M2002 (diamonds), but predicts a flatter 
relation (Tx cx M200 as opposed to Tx oc M^oo)- Lower-mass 
groups have X-ray temperatures in excess of their virial tem- 
peratures (Tvir oc M200 , dotted line), as expected from the 
increase in entropy of these systems. In the highest mass 
group, the presence of cooler reheated gas within the virial 
radius biases its X-ray temperature lower than the virial 
temperature. It is possible that a good fit to higher-mass sys- 
tems will require an even higher heating temperature than 
employed here. 

The dashed and dot-dashed lines in Fig. |3 represent 
power-law fits to the observed temperature-mass relation for 
clusters, compiled by Xu, Jin & Wu (2001), using NFW and 
P model surface brightness profiles respectively. Both fits are 
offset from the simulated relation by up to 40 per cent. A 
likely reason for this offset is that the masses are underesti- 
mated when extrapolating X-ray surface brightness profiles 
to the virial radius. M2002 demonstrated that by estimat- 
ing cluster masses this way (assuming a relation between 
the observable extent of an object and its true size) the re- 
sulting temperature-mass relation was in much closer agree- 
ment with the observations. This conclusion was strength- 
ened when the same dataset was shown to be in good agree- 
ment with the temperature-mass relation for a sample of 
high-resolution clusters observed with the Chandra satellite 
(Allen et al. 2001), for masses measured within a density 
contrast of 2500 (Thomas et al. 2002). 

For groups, the problem of estimating masses from the 
X-ray data is exacerbated by their lower surface brightness, 
and so galaxy velocity dispersion is commonly used as an 
independent measure of the dynamical size of the halo from 
its temperature. We plot the relation between X-ray tem- 
perature and dark matter velocity dispersion, (tdm, for our 
fiducial model in Fig. 2] (Numerical resolution limits the 
number of galaxies per halo to a few in most cases, pre- 
venting us from reliably estimating the galaxy velocity dis- 
persion. Previous simulations have shown that the galaxies 
trace the potential reasonably well, see e.g. Springel et al. 
2001; Dave et al. 2002.) 

The dot-dashed line shows the expected relation if the 
gas contained the same specific energy as the dark matter, 
i.e. Tx oc (Jqi^j. Again, the larger X-ray temperatures of the 
smaller groups flattens this relation, resulting in Tx oc ctdm 
(solid line). To quantify this effect, we plot the ratio of spe- 
cific energies between the hot gas and the dark matter 
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Figure 4. X-ray temperature against dark matter velocity disper- 
sion for groups in the fiducial model (stars) . Square points are for 
when emission is excluded from within the cooling radius (with 
the solid line representing a least-squares fit to these objects). 
Diamonds are from the Radiative simulation studied by M2002. 
The dot-dashed line is the expected relation when the gas and 
dark matter have equal specific energies (T oc 0-2 ), Dotted and 
dashed lines represent best-fits to the group data in Helsdon & 
Ponman (2000) and the compiled group and cluster data in Xue 
& Wu (2000). 
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Figure 5. The ratio between X-ray temperature and dark matter 
velocity dispersion, f3spec, versus virial temperature for groups 
in the fiducial model (stars). Diamonds are from the Radiative 
simulation studied by M2002. 



* The cooling radius is defined as that which encloses hot gas 
with an average cooling time equal to 6 billion years 
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Figure 6. Luminosity-temperature relation for groups in the fidu- 
cial model (stars). Square points are the result when emission 
from within the cooling radius is excluded. Diamonds are from 
the Radiative simulation studied by M2002. Dashed lines repre- 
sent the best-fit relations to the group and cluster samples com- 
piled by Xue & Wu (2000). The dotted line illustrates the best-fit 
power law to the sample of groups found by Helsdon & Ponman 
(2000). 

Pspec — , {I ) 

in Fig. 1^ The smaller objects contain, on average, X-ray 
luminous gas that is nearly twice as hot as the dark matter. 
For the higher-mass objects (fcTvir > IkeV) however, the two 
components are approximately at the same temperature. 

Also plotted in Fig|l|are the best-fits to the group data 
presented in Helsdon & Ponman (2000; dotted line), and the 
group and cluster data compiled by Xue & Wu (2000; dashed 
lines). The Helsdon & Ponman relation is much flatter than 
the Xue & Wu line, due to their sample containing several 
low-temperature groups with low velocity dispersions. Re- 
moving these outliers from the sample would result in closer 
agreement with Xue & Wu's relation. Our simulated rela- 
tion reproduces Xue & Wu's observational determination, 
Tx oc tjgai^ , reasonably well. 

3.3 X-ray luminosity 

The increase in entropy decreases the X-ray luminosity of 
groups. Since this has a bigger effect on less massive objects, 
it leads to a steepening of the — and L^ — a relations, 
which we will quantify here. 

As in M2002 we estimate bolometric X-ray luminosity 
by applying a correction to the soft-band emissivity 

, _ Aboi(rx) ^ mipiAso!t{Ti, Z) 
" " Aeoft(Tx) 2^ itimny ' ^ ' 

i 

where Tx is calculated using equation^] Aboi is the bolomet- 
ric cooling function used in the simulations (see Section 2). 
Figs. |S|&|7| illustrate the Lx — Tx and Lx — a relations 



Figure 7. Lx — ctdm relation for groups in the fiducial model 
(stars). Square points are the result when emission from within 
the cooling radius is excluded. Diamonds are from the Radiative 
simulation studied by M2002. Dashed lines represent the best-fit 
relations to the group and cluster samples compiled by Xue & 
Wu (2000). The dotted line illustrates the best-fit power law to 
the sample of groups found by Helsdon & Ponman (2000). 

respectively, for the fiducial simulation. Again, the model 
agrees reasonably well with both M2002's result and the ob- 
served relations, given the uncertainty in the latter. Further- 
more, increasing the resolution of the smaller objects (which 
are less well-resolved than the larger groups) increases their 
temperature, which will bring the Lx — Tx relation into bet- 
ter agreement with the observed relation (see Section |5J|. 
Note that the "break" in the relation at around 1 keV is 
more prominent in the Lx — Tx relation than in the Lx — cr 
relation, due to the additional effect of entropy increase on 
temperature, causing the relation to steepen further. 

3.4 Baryon fractions 

Although cooling and heating (feedback) are both mecha- 
nisms for increasing the core entropy of groups, they have 
opposing effects on the stellar baryon fraction. Cooling acts 
to increase and feedback decrease galaxy masses, hence 
matching the observed value requires a balance between the 
two processes. Cooling and feedback can also have different 
effects on the hot gas fractions. Although cooling reduces the 
hot gas content, feedback can cause an increase (reheated 
gas that has not moved beyond the virial radius) or a de- 
crease (transfer of energy between outflows and hot halo gas 
increases the fraction expelled from the halo). 

In Fig. |H|we plot mass fractions of baryons, in units of 
the global value {Qb/flm), against the virial temperature of 
each halo for the fiducial model (filled symbols) and M2002's 
Radiative simulation (open symbols) . The total baryon frac- 
tions are plotted in panel (a). Both models appear to be con- 
verging towards a value close to unity at high virial temper- 
atures, but the baryon fractions decrease much more rapidly 
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Figure 8. Baryon fractions, in units of the global value, against 
virial temperature for groups in the fiducial simulation (filled sym- 
bols) and in the Radiative simulation studied by M2002 (open 
symbols). Total baryon fractions are plotted in panel (a), sepa- 
rated into hot gas (diamonds) and cold gas -I- stars (triangles) 
in panel (b). The shaded region represents the 2 global cooled 
fractions measured by Cole et al. (2001) for the Salpeter (upper 
horizontal line) and Kennicutt (lower line) IMF. 

at lower temperatures when feedback is included. The dif- 
ference between baryon fractions reflects the total amount 
of gas expelled from the halo by feedback. 

The baryon fractions are split into hot gas (T > 10"'K, 
diamonds) and cold gas -I- star (triangles) fractions in panel 
(b). The same trend is present between the two models 
for the cooled baryon fractions (since galaxy formation ef- 
ficiency is higher in smaller objects, where average cooling 
times are shorter), but the values at a given temperature 
are different. The cooled fraction is twice as high in the 
M2002 objects, and could be higher: the absence of feed- 



Table 1. Simulations used to study the effects of feedback. From 
left to right, each column gives the model studied, feedback effi- 
ciency, reheating temperature, global fraction of baryons in stars 
at z = and whether the model provides a reasonable match to 
the Lx — Tx relation. 



Model 


e 


fcTsN/keV 


/star 


Lx - Tx? 


Radiative 


N/A 


N/A 


0.15 


YES 


Feedback 


1 


0.1 


0.10 


NO 


Feedback 


0.1 


2 


0.12 


YES 


Feedback 


1 


2 


0.05 


YES 


Feedback 


2 


2 


0.03 


YES 



back in this simulation means that only numerical resolu- 
tion limited these values. Feedback also steepens the rela- 
tion between the hot gas fraction and virial temperature. 
As we discuss in Section 4, objects with virial temperatures 
above the feedback temperature, kTsN, retain a fraction of 
reheated gas which slightly increases their hot gas fraction. 
In smaller systems however, feedback has a larger impact on 
the hot gas, expelling a significant fraction from the halo. 

The shaded region in Fig.|Hlrepresents measured global 
cooled gas fractions (the global value is the asymptotic limit 
for large halo mass), combining the stellar fractions mea- 
sured by Cole et al. using data from the 2MASS and 2dF 
galaxy surveys, with a further 10 per cent of the mass in 
cold gas (Balogh et al. 2001). The upper horizontal line rep- 
resents the value determined when using the Salpeter IMF 
(as used in our feedback model) and the lower for the Kenni- 
cutt IMF. Globally, only just over 3 per cent of the gas forms 
stars in the fiducial model by 2 = 0, lower than the observed 
value by a factor of ~ 3. Increasing the resolution does not 
significantly increase the cooled fraction (Section]^ because 
feedback acts to stabilize this quantity and our simulation 
has already resolved the majority of the galaxy mass density. 
As we discuss in the next section however, a higher cooled 
fraction can be achieved by choosing a lower e which does 
not significantly affect the scaling relations. Additionally, in- 
creasing the metallicity of cooling gas will increase the frac- 
tion of baryons in stars. The metallicity of high-density re- 
gions has been underestimated in our simulation since stars 
in galaxies at the present day have been enriched to solar 
abundance (Edmunds & Phillips 1997), whereas ours have 
Z — 0.3^0. A higher metallicity will increase the fraction of 
stars in a given galaxy. 

Measuring cooled fractions in groups and clusters is dif- 
ficult, mainly due to the lack of sensitivity in the X-ray in- 
struments to measure hot gas masses out to the virial radius. 
Current attempts, however (see for example the compilation 
of data by Balogh et al.), measure values between 10-20 per 
cent at 1 keV. These values are marginally consistent with 
our fiducial result but again, a lower e and/or higher metal- 
licity would improve the match. 



4 A CLOSER LOOK AT THE MODEL 

Given the ability of our fiducial model to adequately repro- 
duce many of the key X-ray observations of groups, we now 
investigate how these results depend on the strength of the 
feedback. To facilitate this task we have run a set of smaller 



The impact of galaxy formation on X-ray groups 9 



ALx/ix = 1 — Lx (feedback) /Lx (radiative). 



(10) 
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Figure 9. Fractional difference in X-ray temperature (stars) and 
luminosity (squares) between matched haloes in the simulations 
with and without feedback, as a function of virial temperature. 
Dotted lines are a guide to the eye. 



simulations (Table with N = 2 x 64^ particles within a 
comoving box of length 25/i~^Mpc, thus keeping the reso- 
lution constant. 

After generating group catalogues for each simulation, 
we also compare the properties of the haloes in two models 
on an object-by-object basis. To do this we first found the 
closest pair, then rejected those which had masses difi^ering 
by more than 15 per cent or separations greater than 20 per 
cent of the average virial radius. These criteria find a good 
match for the majority of the haloes in the two catalogues. 



4.1 Cooling versus heating 

The first comparison we make is between the fiducial model 
(which has strong feedback) and a run with no feedback at 
all (radiative model) . This allows us to quantify the relative 
effects of heating and cooling on the X-ray properties of the 
groups. In particular, we are interested in knowing if the 
model behaves like that proposed by Voit & Bryan (2001), 
i.e. cooling generates the entropy excess by removing hot 
gas, which then either forms stars or is reheated and trans- 
ported to larger radii. If the fiducial model does behave in 
this way, we expect the feedback only to reduce the fraction 
of galactic material, leaving the hot gas fraction (and hence 
X-ray properties) intact. 

Indeed, the Lx — Tx relation for the radiative simulation 
is similar to the fiducial feedback relation, but the former 
is flatter. This effect is quantified in Fig. |5] where we plot 
the fractional difference in X-ray temperature (stars) and 
luminosity (squares) between matched haloes in the two runs 
against their average virial temperature. These differences 
are defined as 



ATx/Tx = 1 - Tx (feedback) /Tx (radiative), 
and 



(9) 



There is no strong variation in the difference between X-ray 
temperatures with system size, although there is a correla- 
tion in the X-ray luminosity difference. For kT^ir ~ O.lkeV 
the objects in the radiative simulation are nearly twice as 
luminous as those in the feedback run, but the opposite is 
the case for systems which are a factor of 10 hotter. (Note 
that because of the small size of the simulations, there is 
only one halo with fcTvir > IkeV.) 

On closer inspection, the trend in luminosity difference 
is down to two effects caused by feedback. The first, which 
particularly affects low-mass systems, is a reduction in the 
mass of hot gas in groups in the feedback simulation which 
was inside the virial radius in the radiative simulation. This 
can be attributed to the increase in pressure from the out- 
flowing reheated material. Secondly, part of the reheated 
gas, which formed stars in the radiative simulation, is still 
within the virial radius in the feedback simulation. As the 
gravitational energy of the halo becomes comparable to that 
given to reheated gas (i.e. Tvir — > Tsn), the second effect 
dominates over the flrst, eventually causing the feedback to 
increase the X-ray luminosity of a group. 

To summarize, we find that our model does behave like 
that proposed by Voit & Bryan (2001) if Tsn ~ Tvir- If 
TsN 3> Tvii then feedback also removes hot halo gas (in 
addition to expelling reheated material) and so reduces the 
group's X-ray luminosity further. Conversely, if Tsn <C Tvir, 
some reheated gas does not gain enough entropy to escape 
and so boosts the group's X-ray luminosity. 



4.2 Varying the model parameters 

We now investigate the effects of varying the feedback model 
parameters, e and Tsn on the X-ray properties of groups. 
Given the simplicity of the model, and the results from the 
previous comparison with a radiative simulation, we should 
be able to predict how varying the parameters affect the 
results. 

First of all, we consider the effect of varying the feed- 
back parameters on star formation efficiency. The star for- 
mation rate, which determines the amount of energy and 
metals supplied to the interstellar medium, varies inversely 
with /heat (equation |3Jl, since more reheated cold gas leaves 
less material available to form stars. Note /heat oc e/fcTsN, 
so we would expect a smaller energy budget and/or a higher 
heating temperature to increase the star formation rate. 

Panel (a) of Fig. 1101 illustrates star fractions against 
virial temperature for groups in the radiative model plus 
feedback models with varying parameter values (given in 
the legend). As expected, the star fraction increases when 
decreasing e for fixed fcTsN, and the model with e = 0.1 con- 
tains approximately the correct amount of cooled material 
(~ 10 per cent, c.f. Fig. |HJ. The impact of varying e on the 
total baryon fraction is shown in panel (b). As e is increased, 
more gas is reheated (at the expense of forming stars) and 
expelled from each halo, decreasing its baryon fraction. 

Contrary to expectations however, the star fraction in- 
creases when fcTsN is decreased from 2 to 0.1 keV, with 
the total baryon fraction of the low temperature case being 
similar, on average, to the radiative model. As was found for 
K2002's thermal feedback model, if Tsn < Tvir for the haloes 
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Figure 10. Fraction of baryons in stars (panel a) and total baryon 
fraction (panel b) against halo virial temperature for the radiative 
run and the runs with varying feedback parameters. 



being studied (as is the case for the kTsN ~ 0.1 keV model), 
a substantial fraction of reheated gas is not hot enough ^ 
to escape from the gravitational potential, and will instead 
re-cool and boost the star formation rate. 

The effect of feedback in the fiducial model was to 
slightly steepen the Lx — Tx relation that was produced by 
radiative cooling alone. Fig. 1111 illustrates the effect of vary- 
ing the parameters on the Lx — Tx relation. Although the 
simulations do not contain enough objects to perform a de- 
tailed study of the shape of the Lx — Tx relation, we see 
from our present results that the overall normalization of 
the relation is fairly insensitive to e when Tsn > Tvir, except 



Note however that some gas will have been hot enough to escape 
from progenitor haloes at earlier times 



Figure 11. Lx — Tx relations for the radiative simulation and 
runs with varying feedback parameters. 



Table 2. Simulations used to study the effects of varying mass 
resolution 



Model 


N 




€p/h ^kpc /cool 


LR 


2 X 32^ 


5.2 X 10^ 


10 0.031 


HR 


2 X 64^ 


6.5 X 10* 


5 0.038 



that the relation becomes slightly steeper as e is increased 
from (radiative) to 2. Although the heating temperature is 
the same in these runs, increasing e increases the amount of 
outflowing material, thus having a larger impact on the hot 
gas (as discussed in the previous section). When Tsn < Tvir 
however (as is the case with fcTsN = O.lkeV) the luminosi- 
ties are boosted and temperatures decreased considerably. 
As was observed from the previous figure, reheated mate- 
rial in the latter run is too cold to completely escape from 
a group-sized halo: the lack of buoyancy causes this gas to 
boost the central density, hence boosting its X-ray luminos- 
ity and lowering its X-ray temperature. 

Although the fiducial run provides the best fit to the 
observations (at least for low-mass groups), we will show in 
the next section that the lack of numerical resolution in the 
smaller objects artificially fiattens the relation. 



5 RESOLUTION EFFECTS 

The results presented in this paper have all been for fixed 
resolution. In this section, we present results from a limited 
study of how varying the resolution affects the X-ray prop- 
erties of the groups. Table|5|lists properties of the two sim- 
ulations studied here, a low-resolution simulation (labelled 
LR) with the same resolution used for the runs studied pre- 
viously, and a high-resolution simulation (HR) with 8 times 
as many particles. Both runs were performed for a comov- 
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Figure 12. Lx — Tx relation for the groups in run LR (squares) 
and HR (stars). Matched objects are connected with solid lines. 
The N = 2x 128^ fiducial relation (same resolution as LR) is also 
plotted (diamonds). Numbers beside each point give the number 
of hot gas particles within the virial radius of the object. Solid 
and dashed lines illustrate best-fit observed relations presented 
by Helsdon & Ponman (2000) and Xue & Wu (2000) respectively. 



ing box-size of 12.5 h~^Mpc and contain the same large-scale 
power. Although such a box is too small for quantitative pre- 
dictions, it allows us to compare a small sample of the same 
objects at different resolution, without requiring a large in- 
vestment in CPU time. The fiducial feedback parameters 
were adopted, namely e = 2 and fcTsN = 2keV. 

The final column in Table|5|lists the cooled baryon frac- 
tion for the 2 runs; a factor of 8 increase in resolution only 
increases this fraction by 20 per cent. This is because 
the majority of the mass in galaxies is already resolved at 
the lowest resolution (i.e. 32 or more particles, the lowest 
number than can cool efficiently.). Furthermore, the feed- 
back acts to stabilize gas cooling since the heating rate is 
directly proportional to the cooling (star formation) rate. 
As was stated previously, an overall increase in the cooled 
fraction must therefore come from choosing a lower value of 
e and/or higher metallicity of cooling gas. 

All runs resolve 3 common objects; the corresponding 
Lx — Tx "relation" for each simulation is plotted in Fig. 1121 
Particles from within the cooling radius of each object are 
excluded from the calculation (which affects only the lumi- 
nosity of the largest object). Note that the luminosities are 
more or less unchanged for all 3 objects when the resolution 
is varied, due to the stabilizing effect of feedback. For the 
largest object, which contains between 2800 (LR) and 18,000 
(HR) hot gas particles, its X-ray temperature increases only 
a little, unlike the smaller 2 objects (between 100 and 1000 
particles), where it increases significantly. We therefore ex- 
pect the slope of the Lx — Lx relation to steepen when the 
emission-weighted temperatures for all objects have con- 
verged. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we examined the influence of the galaxy pop- 
ulation on the X-ray properties of groups at z = that form 
in a cosmological simulation. In particular, we were inter- 
ested in determining how including feedback of energy from 
supernovae in the simulations (which decreases the star for- 
mation rate) affects the X-ray properties of the groups. Mo- 
tivated by the "break" in the X-ray luminosity-temperature 
relation at the interface between group and cluster scales 
(fcTx ~ 1 — 2keV), we devised a feedback model based on 
heating the gas to a fixed temperature, independent of the 
available energy budget. A given heating temperature, kTsN, 
and energy budget, e, then determines the star formation 
rate and the corresponding mass of reheated gas. 

For our main results, we studied a model with strong 
feedback (e = 2, twice the available energy from supernovae) 
and a heating temperature, fcTsN = 2keV, comparable to 
the Lx — Lx break. This model generated the required ex- 
cess entropy in the cores of the groups and so was able to 
adequately reproduce the X-ray scaling relations over the 
range of temperatures appropriate to groups {kT < 2keV). 
However, the fraction of baryons in galaxies was around a 
factor of 3 lower than observed (assuming a Salpeter IMF). 

We then compared our fiducial feedback model to a 
model with no feedback, and found that the Lx — Tx re- 
lations were similar in normalization, but the feedback re- 
lation was slightly steeper. Radiative cooling acts to reduce 
the hot gas fraction (and hence X-ray luminosity) in groups, 
removing low-entropy material (which forms stars) and caus- 
ing higher-entropy gas to flow in to replace it. In low-mass 
groups (where TVir < Lsn), feedback reduces the luminosity 
further because the outflowing reheated material increases 
the pressure of the hot gas. However, as Tvir Lsn, the ef- 
fect of feedback on the hot gas diminishes due to the group's 
increasing ability to retain a larger fraction of material. As 
a result, some of the gas that was turned into stars in the 
radiative run is still within the virial radius in the feedback 
run, leading to an increase in the luminosities of high-mass 
groups. 

The fraction of reheated material per mass of stars 
scales as e/fcTsN, and so the star formation rate ought to 
vary inversely with this parameter combination, since more 
reheated gas leaves less material available to form stars. This 
is only true if the heating temperature is comparable to or 
larger than the virial temperature of the systems being stud- 
ied. For such a case, decreasing e allows a reasonable match 
to both the global fraction of baryons in stars and the X- 
ray properties of the groups. However, when the heating 
temperature is low compared to the virial temperature of 
groups (we studied the case when fcTsN = O.lkeV), reheated 
gas cannot escape from the gravitational potential and in- 
stead re-cools and forms more stars. As a result, this excess 
of cooling gas decreases the core entropy of the groups and 
consequently boosts their luminosity. 

Finally, increasing the resolution of our simulations ap- 
pears to have little effect on the X-ray luminosities of groups 
but the temperatures of the smaller, poorly-resolved objects 
increased. Based on these results, numerical convergence for 
the range of group masses will act to steepen the Lx — Tx 
relation. 
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